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ABSTRACT 

In 1996, a major radio flux-density outburst occured in the broad-line radio galaxy 3C 111. It was followed 
by a particularly bright plasma ejection associated with a superluminal jet component, which has shaped the 
parsec-scale structure of 3C 111 for almost a decade. Here, we present results from 18 epochs of Very Long 
Baseline Array (VLBA) observations conducted since 1995 as part of the VLBA 2 cm Survey and MOJAVE 
monitoring programs. This major event allows us to study a variety of processes associated with outbursts of 
radio-loud AGN in much greater detail than has been possible in other cases: the primary perturbation gives rise 
to the formation of a leading and a following component, which are interpreted as a forward and a backward- 
shock. Both components evolve in characteristically different ways and allow us to draw conclusions about 
the work flow of jet-production events; the expansion, acceleration and recollimation of the ejected jet plasma 
in an environment with steep pressure and density gradients are revealed; trailing components are formed in 
the wake of the primary perturbation possibly as a result of coupling to Kelvin-Helmholtz instability pinching 
modes from the interaction of the jet with the external medium. The interaction of the jet with its ambient 
medium is further described by the linear-polarization signature of jet components traveling along the jet and 
passing a region of steep pressure/density gradients. 

Subject headings: galaxies: individual: 3C1 1 1 - galaxies: active - galaxies: jets - galaxies: nuclei 



1. INTRODUCTION 

Direct evidence for the existence of bulk relativistic out- 
flows along the jets in blazars and other radio-loud ac- 
tive galactic nuclei (AGN) comes from Very-Long-Baseline 
Interferometry (VLBI) observations. The first evidence 
for apparently superluminal structural changes was found 
from ch anges in the fringe visibility curves o f 3C279 and 
3C273 (IWhitnev et alJ[l97llCohen et al]fl97lh . Subsequent 
higher-quality VLB! obse rvations (see, e.g., compilation by 
Verm eulen & Cohenlll994i and references therein) have es- 
tablished the "core-jet" type milliarcsecond-scale structure 
of compact extragalactic jets: the core being a bright and 
unresolved flat-spectrum component at the end of a lin- 
ear structure, and the jet being composed out of individ- 
ual steep-spectrum components or "knots". The knots fre- 
quently move away from the core with apparent veloci- 
ties exceeding the speed of li ght. Monitoring observa- 
tions of large source samples jV ermeulen & Cohen Il994t 
Uorstad et al.1l200UlKetlermann et aLll2004t iPiner et al.ll2007h 
have provided important statistical tools for probing relativis- 
tic beami ng and the intrins ic properties of extragalactic ra- 
dio jets (ICohen et al.1 120071) . their intrinsic brightness tem- 
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peratures ( Ho man et al J 120061). or their Lorentz factor dis- 
tribution ( Kellerma nn et al mool and luminosity function 
(ICara & Listerll2007l) . 

The relativistic-jet model (e.g., Blan dford & Konigll 19791) 
has become the de-facto paradigm in multiwavelength re- 
search on blazars and other AGN, but VLBI observations 
have demonstrated that the basic concept of ballistically- 
moving isol ated jet knots is clearly ove rsimplified: jet cur- 
vature (e.g., Vermeulen & Cohen 1994), stationary compo- 
nents (e.g.. iJorstad et al.ll2001l). and non-r adial and acceler- 
ated motions (e.g.. iKellermann et alj|2004l) . are found to be 
common features of relativistic jets. Within individual jets, 
there are often characteristic velocities suggesting the pres- 
ence of an underlying continuous jet flow, but the "com- 
ponents" themselves most likely represent patterns moving 
at a different speed than the under lying flow, e.g., as hy- 
drodynamically pro pagating shocks dMarscher & Gear|[l985l : 
iHughes et alJ[l985l) . 

Recent years have brought major impro vements in nu meri- 
cal simulations of relativistic jets (see, e.g., [Gomez 2005], for a 
review). It is now possible to simul ate three-dimensional rel- 
ativistic jets (e.g.. lAlov et al.ll2003l) and t o compute the rela- 
tivistic processes (e.g., iGomezetal] 19971) that transfer hydro- 
dynamic results into observed brightness distributions (e.g., 
relativistic light abberation and light travel time delays). In 
particular, interactions between strong perturbations or shocks 
with the unde rlying jet flow an d the jet-ambient medium can 
be simulated (Agudo et al. 2001). With these new techniques, 
it is now possible to compare the generation, propagation and 
evolution of emission features in simulated and observed rel- 
ativistic jets. 

The nearby (z=0.049) 9 broad-line radio galaxy 3C 111 
(PKS B 0415+379) shows a classical FRII morphology on 
kiloparsec-scales spanning more than 200" with a highly col- 
limated jet connecting the central core and the northeastern 

9 Assuming H = 7110ns" 1 Mpc~', fi M = 0.3, Q A =0J (lmas = l.Opc). 
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lobe in position angle 63° whi le no counter jet is obser ved to- 
wards the southwestern lobe (Lin field & Perleylll984t) . This 
asymmetry is usually explained via relativistic boosting of 
the jet and de-boosting of the counter-jet. 3C 111 exhibits 
the brightest compact radio core at cm/mm wavelengths of 
all FRI I radio galaxies, a blazar-like spectral energy distri- 
bution (Sgueraet al. 2005), and it was one of the first (and 
only) r adio galaxies in which superluminal motion was de - 
tected (lGoetzetal.ll 19871: iPreuss. Alef. & KeUermannlll988l) . 
Moreover, the (sub-) parsec scale j et of 3C 11 1 is in timately 
related to its high-energy emission: lMarscherld2006l) reports 
a disk-j et connection, similar to the well-established one in 
3C 120 dMarscher et al.ll2002l) . in the sense that dips in the X- 
ray light curve indicate accretion events which are followed 
by VLBI jet component ejections. Recently, R.C. Hartman 
& M. Kadler (in prep.) showed that the gamma-ray source 
3EGJ04 16+3650 can be decomposed into multiple individ- 
ual sources inside the EGRET full-band point-spread func- 
tion, revealing a significant signal from the nominal position 
of 3C 1 1 1 in the higher-resolution, high-energy band above 
1 GeV. This association of 3C 1 1 1 w ith 3EGJ04 16+3650 
whi ch had origina ll y been suggested bv lHartman et al.l (1 19991) 
and ISguera et al.l d2005), makes 3C 111 one of the very 
rare radio galaxies detected at gamma-ray energies and sup- 
ports the view that this source may be considered a lower- 
luminosity version of powerful radio-loud quasars. 

Here, we report the results from ten years of Very-Long- 
Baseline Interferometry (VLBI) observat ions of 3C 1 1 1 
part of the VLB A 2 cm S urvey 10 (iKellermann et al. 
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JAVE" dLister & Homanll2005t iHoman & Lister! 120061) . We 
investigate the parsec-scale source structure during a major 
flux-density outburst and during its aftermath. We find 
that this outburst was associated with the formation of an 
exceptionally bright feature in the jet of 3C 1 1 1 . A variety of 
processes (beyond the predictions of simple ballistic motion 
models) are observed and discussed in view of modern 
relativistic-jet simulations. In Sect. [2] our observations and 
the data reduction are described. A detailed report of the 
observational results is given in Sect. [3] In Sect. @] we discuss 
the various processes observed in the jet of 3C 1 1 1 as a result 
of the outburst and during the propagation of the new jet 
feature along the jet. In Sect. [5] we put these results into the 
context of future simulations and observations with the goal 
of understanding the production mechanisms of AGN jets. 

2. OBSERVATIONS AND DATA ANALYSIS 

3C 111 has been monitored as part of the VLBA 2 cm Sur- 
vey program since April 1995. Th e observational details 
are given by IKellermann et all (119981) . Following the meth- 
ods described there, the data from 17 epochs of VLBA ob- 
servations of 3C111 between 1995 and 2005 (see Table [TJ 
were phase and amplitude self calibrated and the brightness 
distribution was determined via hybrid mapping. An addi- 
tional epoch from June 2000 was made available to us by 
G. Taylor. The polarization calibra tion was performed as de- 
scribed in iLister & Homanl (120051) . Two-dimensional Gaus- 
sian components were fitted in the («,v)-domain to the fully 
calibrate d visibility data of each epoch using the program 
DIFMAP (Shepherd 1997). The parameters of each model fit 
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FIG. 1. — University of Michigan Radio Astronomy Observatory light 
curves of 3C 1 1 1 at 4.8 GHz, 8 GHz, and 14.5 GHz. The shaded areas indicate 
the ejec tion epochs of the individually labeled jet components as discussed 
in Sect. 13.21 The lightest shading corresponds to minor ejections of the rel- 
atively weak components I,J,K,L with flux densities S below 0.2 Jy, medium 
shading corresponds to components B,C,G,H,M,N with 0.2Jy < S < 0.6Jy 
and the darkest shading to components E and F with S > 0.6Jy. 
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FIG. 2. — Spectral-index curves of 3C 111 between 14.5 GHz and 8 GHz 
(top), and 8 GHz and 4.8 GHz (bottom) from the UMRAO monitoring pro- 
gram. The shaded areas are the same as in Fig.fT] 

at the various epochs are given in Table [2] The models were 
aligned by assuming the westernmost component (namely, the 
"core") to be stationary so that the position of jet components 
can be measured relative to it. Because of the coupling of 
the flux densities of nearby model components, the uncertain- 
ties in the component flux densities are larger than the formal 
(statistical) errors unless the given model component is far 
enough from its closest neighbor. Throughout this paper, er- 
rors of 15 % are assumed for the flux densities of individual 
model-fit components. In most cases, this should be consid- 
ered a conservative estimate that accounts for absolute cal- 
ibration uncertainties and fo rmal model-fitting uncertainties 
(see, e.g.. lHoman et al.l[2002l) . Position uncertainties were de- 
termined internally from the deviations of the data from linear 
motion. 

3. RESULTS 
3.1. The 1996 Radio Outburst of3Clll 
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FIG . 3 . — Naturally weighted images of the parsec-scale jet of 3C 1 1 1 from the 2 cm VLBA monitoring. A common restoring beam of (0.5 X 1 .0) mas at 
P.A. 0° was used. The total recovered flux density in each image, the rms noise, and the lowest contours for each image are given in Tableff] Contours increase 
logarithmically by a factor of 2. Only components E, G, H and their corresponding trailing components are indicated by circles enclosing a cross. 



A strong flux density outburst of 3C 1 1 1 occurred in 1996, 
which was first visible in the mm band and some months 
later at lower radio frequencies. This outburst was first de- 
tected at 90 GHz with the IRAM interferometer at Plateau 
de Bu re in January 19 96 with flux densities greater than 
10 Jy (lAlefetal.lfT998h . at 37 GHz in March 1996, and at 
22 GHz in August 199 6 with the Metsahovi radio observatory 
(Terasranta et al. 2004). Figure [T] shows the single-dish radio 
light curves of 3C 1 1 1 at 4.8 GHz, 8 GHz, and 14.5 GHz ob- 
tained from the UMRAO radio-flux -density monitoring pro- 
gram ( Aller. Alien & Hughes! 120031) . These data show that 
from early 1996 on, the radio-flux density of 3C111 was 
rising at 14.5 GHz, reaching its maximum in late 1996. At 



the two lower frequencies, the flux-density maximum was 
reached at subsequent later times, in mid 1997 at 8 GHz and 
in late 1997 at 4.8 GHz. The profile of the outburst in the 
flux-density vs. time domain shows a narrow, high-amplitude 
peak between early 1996 and late 1997, which is almost sym- 
metric. After late 1997, a slower-decreasing component dom- 
inates the light curves, most clearly visible at 14.5 GHz. 

Figure [TJ shows that the flare propagated through the 
spectr um as qualitatively ex pected by standard jet theory 
(e.g jMarscher & Gearll985l) : high-frequency radio emission 
comes from the most compact regions of the jet, the emission 
peak shifts to lower frequencies as a newly ejected jet com- 
ponent travels down the jet and becomes optically thin. The 
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peak flux density shifted with frequency at about 10 GHz yr 1 . 

The evolution of the spectral index, a {S ~ v a ), for 
(14.5/8.0)GHz and (8.0/4.8)GHz is shown in Fig. El Before 
1996, the sampling was too sparse to derive the change of the 
spectral index in the (14.5/8.0) GHz band. Between, 8.0 GHz 
and 4.8 GHz, the spectral index was approximately -0.7 dur- 
ing the pre- 1996 period. The radio flux-density outburst in 
1996 corresponded to a subsequent flattening of the spectrum 
with a maximum spectral index, a ~ 0, in the (14.5/8. 0)GHz 
band reached in mid 1996. In the post-outburst period be- 
tween 1998 and 2004, a was typically in the range -0.5 
to -0.7 between 14.5 GHz and 8.0GHz and slightly steeper 
(-0.7 to -0.9) in the (8.0/4.8) GHz band. The overall steeper 
spectral index at lower frequencies can be understood as the 
contribution of optically thin large-scale emission from the 
radio lobes of 3C 1 1 1 to these single-dish light curves. 

3.2. VLBA Monitoring Results 

Figure [3] shows the variable parsec-scale structure of 
3C 111 at 18 different epochs of VLBA observations between 
1995.26 and 2005.73. The variable source structure can be 
described by a classical one-sided core-jet morphology in the 
first two epochs with typical velocities of the outward moving 
jet components of about 1.4 to 1.7 mas yr" 1 corresponding to 
about 5 c. In 1996.82 a new jet component, even brighter than 
the core, dominated the source structure. By 1997.19, this 
new component was even brighter (~ 3.4Jy) and in the fol- 
lowing epochs it traveled along the jet while it became gradu- 
ally more stretched out along the jet-ridge line. 

Model Fitting: — In Fig.@] the radial distance of the various 
model fit components from the core is shown as a function of 
time. The component identification was based on a compari- 
son of the positions and flux densities, and a linear regression 
of the distances from the core as a function of time was used 
to determine the kinematics. The derived component veloci- 
ties are tabulated in Table [3] The early outer jet components 
(A, B, C, D) of the 1995.26 epoch can be traced over two to 
four epochs before their flux densities fall below the detection 
threshold (compare Fig.|5}. In late 1996 and early 1997, the 
source structure was dominated by the emission of the core 
and the newly formed jet components E and F, with E being 
the leading component. The two components traveled out- 
wards with a mean apparent velocity of (1.00 ±0.02) mas yr" 1 
and (0.64 ± 0.07) mas yr" 1 , respectively. Before mid 1997, 
component F was substantially brighter than component E 
but after that, its flux density dropped steeply. F was not de- 
tected at any epoch later than 1998. 18, while E was still about 
800 mJy at that time. The light curves of E and F reproduce 
qualitatively the two-component shape of the flux-density out- 
burst in Fig. Q] with component F being responsible for the 
narrower and higher-amplitude peak between early 1996 and 
late 1997 and component E dominating the slower-decreasing 
tail of the outburst after late 1997 (compare Fig. [5] and dis- 
cussion below). In the following epochs, E split into four dis- 
tinct components (E 1, E2, E3, E4) at distances of 3.5 mas 
to 4.5 mas from the core. E2, E3, and E4 all moved at sub- 
sequently slower speeds than E 1, resulting in an elongated 
morphological structure of the associated emission complex. 

In later epochs, new components have been ejected from 
the core into the jet. The two strongest components (G,H) 
can be traced through the following eight and nine monitoring 
epochs, respectively. Component H split into three individual 
components in 2004.27 and a fourth associated component 
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FIG. 4. — Core separation of model-fit components vs. time. Crosses 
represent components which could not be cross identified over the epochs. 
Position uncertainties have been estimated from the internal deviations of the 
data from linear motions for each component. The uncertainties for the posi- 
tion of E 1 have been determined separately for the pre-2004 and post-2004 
epochs because of the partial resolution of this component after 2004. Uncer- 
tainties for components with less than three epochs were estimated from other 
components at similar positions in the jet and with similar flux densities. 




FIG. 5. — Flux-density evolution of the core and the jet components with 
time. For clarity, only components ejected before 2001.50 are shown. The 
flux densities of E 1, E2, E 3, and E4 were added for the post- 1999 epochs 
and a flux-density weighted effective position was calculated to display the 
flux-density evolution of the blended feature that would be visible at lower 
resolution. Note that components E4 and G are blended in epoch 2004.27 
and that the flux density of E 4 may be overestimated for this epoch. 

was seen from 2004.80 on. In the following, we refer to the 
components E 1, and H 1 as the "leading components" and to 
E2, E3, E4, H2, H3, and H4 as the "trailing components" 
of E and H, respectively. 

For the pc-scale jet of 3C 111, the ejection epochs of the 
individual jet components can be determined from the linear 
regression by back-extrapolating the component trajectories 
to the core. In Fig. Q] these ejection epochs and the associ- 
ated uncertainties are indicated as shaded areas. It is apparent 
that the ejection of the components E and F coincides with 
the onset of the major flux-density outburst in 1996 described 
above. The following major component ejections (G, H, and 
the combined M/N event) all have direct counterparts in local 
maxima of the radio light curve, especially at 14.5 GHz. Fig- 
ure|2]shows that all these ejection epochs coincided with local 
maxima of the spectral index in the 14.5/8.0 GHz band. Be- 
tween 2002 and 2004, a number of minor component ejections 
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FIG. 6. — Flux-density evolution of component E and its trailing compo- 
nents vs. the distance traveled from the core. Note that component E4 is 
blended with component G in epoch 2004.27 and that its flux density may be 
overestimated for this epoch. 
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FIG. 7. — Flux-density evolution of component G and H and its trailing 
components vs. distance traveled from the core. 



took place but the regression-fit quality (due to the neareby 
components, the low flux densities and the small time base- 
line) only moderatly constrains the ejection epochs. In addi- 
tion, the time sampling of UMRAO observations in this time 
range is relatively poor, in particular from mid 2001 to mid 
2003. 

Flux Density Evolution: — Figure[5]shows the brightness evo- 
lution of the core and the jet components that have been 
ejected prior to 2001.5. Apparently, the trailing components 
E4 and H4 appeared first in a rising state, i.e., they first in- 
creased in flux density before they became fainter in later 
epochs. Component F showed an extraordinary steep de- 
crease in brightness in 1997-1998. In Fig.|6]and Fig. [7] the 
flux-density evolution of the components E, G, and H and the 
associated leading and trailing components are shown with 
distance traveled from the core, respectively. The ejecta first 
rose in flux density within the inner 1 mas from the core, then 
they showed a decline about almost three orders of magni- 
tude in the following decade, exhibiting a plateau or broad lo- 
cal maximum in 1998-2000 at a distance from 2-4 mas from 
the core. Component H and its leading and trailing compo- 
nents exhibited a similar behavior although on about an order 
of magnitude lower flux-density levels and at slightly further 
downstream, 4-6 mas from the core. Component G, in spite of 
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FIG. 8. — Brightness temperatures of model-fit components as a function 
of their distance to the core. The brightness temperatures of components 
belonging to the E-, F- and H-components are indicated by filled black cir- 
cles. The solid line represents a least-squares fit to all but the E-, F- and 
H-components. The slope of the regression curve is -2.4 ± 0.2. 
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FIG. 9. — Brightness temperatures of component E, its leading and trailing 
components, and component F as a function of their distance to the core. 
The two regimes of brightness-temperature gradient discussed in the text are 
indicated with dashed lines. 

the fact that it does not appear to have split into leading and 
trailing components like E and H, did exhibit a pronounced 
flux density maximum after 2002, as well, about 6 mas from 
the core. 



The 7b gradi e nt alo ng the jet: — Following iKadler et alj 
(2004); IKadler! (120051) . the power-law index s, which de- 
scribes the brightness temperature gradient via 7b oc r\ can 
be parametrized as 



s = l + n + b(\ -a) 



(1) 



where I, n and b are the power law indices that describe the 
gradients of jet diameter d oc r l , particle density n e oc r n , and 
magnetic field B oc r b with distance r from the core, respec- 
tively. Therefore, measuring the brightness temperature gra- 
dient provides a method to constrain the critical physical prop- 
erties along the jet and abrupt changes in the 7t,-gradient can 
highlight regions in the jet where the density, magnetic field, 
or jet diameter change rapidly. 

Figure [8] shows the brightness temperatures of all jet com- 
ponents in the parsec-scale jet of 3C 1 1 1 at 2 cm wavelength 
between 1995 and 2005 as a function of their core distance. 
In general, the brightness temperature of all components de- 
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creased as the components traveled outwards but an approx- 
imation with a simple power law does not yield a good fit to 
the full data set (Xred = l-8> 115 degrees of freedom [d.o.f.]). 
Visual inspection of Fig. [8] shows that this is due to the E-, 
F- and H- components and their leading and trailing compo- 
nents, respectively. This behavior is different than expected 
for a straight and stable jet geometry in which the power- 
law dependences of the particle density, the magnetic field 
strength and the jet diameter on the core distance predicts that 
the brightness temperature along the jet can be described with 
a well-defined power-law index s. Most extragalactic parsec- 
scale jets which do not show pronounced curvature, show a 
power-law decrease with increasing distance from the cor e 
and power-law indices typically around -2.5 (Kadler 2005). 
In fact, excluding the E-, F- and H-components from the fit 
yields a statistically better result (xjL = 1.3, 5 2 d.o .f.) and 
a gradient of -2.4 ± 0.2. In Sect, Iljjand Sect. 03] we dis- 
cuss possible physical reasons for the different behavior and 
nature of these components. The measured relation between 
component sizes and distance along the jet is affected by a 
large degree of scatter and does not provide independent in- 
formation from the flux density and brightness temperature 
plots described above. Therefore we do not show plots of 
component size versus jet distance. 

The brightness-temperature gradient of component E was 
first flat or inverted immediately after the creation of this new 
component within approximately 1 mas from the core and 
then reached steep values of -2.5 to -2.8 (regime I; com- 
pare Fig. |9]l through 1997 when the component traveled from 
1 mas to 2 mas. Between 2 mas and 4 mas, the determination 
of the brightness-temperature gradient requires an identifica- 
tion of component E with either component El or E 3 (see 
below). Independently of this identification, the brightness- 
temperature gradient eventually changed to very steep values 
(s < -5) beyond 5 mas from the core (regime II). Component 
F began its very rapid decline in brightness temperature at 
a very small distance from the core (< 0.7 mas) with an ex- 
tremely steep 7b-gradient (s < -8). 

Linear Polarization: — From 1995 to 2002, 2cm-Survey ob- 
servations were done in left circular polarization only, so 
no linear-polarization information can be derived from these 
data. MOJAVE observations (from 2002 on) are done in full- 
polarimetric mode. Figure [10] shows our polarization data 
through September 2005. 

In October of 2002, component H was about 2.5 mas from 
the base of the jet and showed a fractional polarization of 5% 
to 10% increasing towards the downstream side of the compo- 
nent. The electric vector position angle (EVPA) displayed by 
the component was approximately aligned with the jet. In this 
epoch, the jet material just downstream of component H at 
~ 3.3 mas from the base of the jet was more highly polarized, 
exceeding 20% fractional polarization on the jet's southern 
side, and the EVPA of the polarization turned to be about 45° 
to the main jet direction. 

By August of 2003, component H had entered a region ap- 
proximately 3.3 mas from the core and its observed polariza- 
tion was now similar to the emission in this same region ob- 
served in the previous epoch. The observed fractional polar- 
ization of H now climbed sharply to values in excess of 20% 
toward the jet's southern side while there was no detectable 
polarization from the northern side of H. The observed EVPA 
of H had rotated further to be approximately 60° to the main 
jet direction. However, the EVPA on the southern-most side 
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FIG. 10. — Naturally weighted images of the linear-polarization structure 
of 3C 111 between 2002 and 2005. The restoring-beam dimensions and ori- 
entations for each epoch are indicated by a cross to the left of each Stokes-I 
image. Stokes I contours start at 1 mjy/beam and increase by factors of 2. 
Fractional polarization is over-plotted on the Stokes-I contours in color. To 
the right of the Stokes-I images are the polarization intensity contours starting 
at 1 mJy/beam and increasing by factors v2. The polarization contours are 
over-plotted with tick-marks representing the electric vector position angle. 
A single Stokes-I 1 mjy contour surrounds the polarization image to show 
registration. The dotted line marks the distance of 3.3 mas from the core 
where the most pronounced changes of the polarization properties take place 
(see text). 

of H was approximately perpendicular to the jet. 

After component H passed through this region (epochs 
April 2004 through September 2005), it split into a number of 
subcomponents as described earlier, and its polarization grad- 
ually became more uniform. Consistent fractional polariza- 
tion of 5% to 10% was approached with the electric vectors 
approximately perpendicular to the local jet direction. 

The much weaker component, I, developed polarization 
very similar to H as it passed through the same region, about 
3.3 mas from the core, with fractional polarization exceeding 
20% toward the southern side of the jet with an EVPA at ap- 
proximately 45° to the main jet axis. This is also the same 
region of the jet in which component E had broken up into a 
number of sub-components. In future epochs, we will have 
the opportunity to follow component K as it passes through 
this same region. 

4. DISCUSSION 

In this section, we discuss the aftermath of the major out- 
burst in 3C 111 in 1996 and the following component ejec- 
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tions through 2005. We organize the subsections of our dis- 
cussion according to the downstream distance from the VLBI 
core where we observe the effect of interest. 

4. 1 . Within 1 pc: Forward and Reverse Structures 

Numerical simulations (lAlov et al.l 120031) show that an 
abrupt perturbation of the fluid density at the jet injection 
point during a short time propagates downstream, evolves 
spreading asymmetrically along the jet and finally splits into 
two distinct regions. Both of these two regions have enhanced 
energy density with respect to the underlying jet, and they 
emit synchrotron radiation. The leading (forward shock) and 
the following region (reverse shock) have higher and lower 
Lorentz factors, respectively, than the underlying jet. Thus, 
they should separate with time as they propagate downstream 
in the jet. 

Component F matches the description of a backward mov- 
ing wave associated with the major injection into the jet of 
3C111 after the flux-density outburst of 1996. It follows 
the trail of component E but at a lower speed. If compo- 
nent F is identified with a reverse shock and component E 
with a forward shock, it is possible then to compute the size 
of the shocked region dPerucho et alj|2007b . In 1996.82 and 
1997. 19, E and F were both very bright and separated by only 
~ 0.3 pc in projected distance. During these two epochs, F 
was 300 mJy to 5 00 mJy brighter th an the leading compo- 
nent E. Following lAlov et alT ( l2003l) . a backward shock can 
be brighter than a forward shock if the latter is beamed in a 
cone smaller than the viewing angle due to its larger speed. 
We have examined the Doppler factors of components E and 
F for the range of possible viewing angles (see Appendix lAl 
and the measured velocities and conclude that this alone can- 
not explain the brightness difference between component E 
and F bec ause the differenc e in apparent speed is not large 
enough. Jorsta dTt al.l (120051) point out that backward shocks 
can be brighter than forward shocks as long as the distur- 
bance is prolonged and there is a continuous supply of par- 
ticles entering from the underlying jet through the shock re- 
gion. Within half a year, between 1997. 19 and 1997.66, F lost 
about half of its brightness. This extraordinarily fast dimming 
of the backward shock can be caused by the lack of input of 
particles from behind, i.e., a lower plasma ejection rate after 
the primary inj ection possibly due to a depletion of th e inner 
accretion disk dMarscher et alj|2002 l ; iMarscherl 12006) which 
feeds the plasma injection. 

Component F can also be interpreted as a rarefaction propa- 
gating backwards in the reference frame of the ejected blob of 
gas. A rarefaction is produced when the blob is overpressured 
with respect to the jet, as this overpressure causes the front 
to accelerate in the jet, thus leaving a rarefied region between 
the head of the blob (forward shock) and its rear part, which is 
still slower (it moves with the injection velocity). In this case, 
the emission in component F could be associated to the denser 
and overpressured gas in the blob which has still not been rar- 
efied. This gas would cease to emit as soon as it reaches the 
rarefied region, which may also explain the sudden decrease 
in brightness of this component. An extended discussion on 
the nature of co mponent F and the e volution of its brightness 
will be given in lPerucho et al.l (120071) . 

4.2. Between 2pc and Ape: Expansion and Acceleration 

It is not a-priori clear with which post-split-up component 
the original feature E should be identified after 1999. A natu- 
ral identification would be the leading component El but that 



requires an acceleration of this component (see Fig. 4) from 
/3a PP ,E = 3. 26 ±0.07 to /3 ap p,Ei = 5.5 ±0.1 between 1998.18 and 
1999.38. This may be interpreted in terms of an expansion of 
the jet in a rarefied medium. Taking an angle to the line of 
sight of 19° (see Appendix lAl). the component would be ac- 
celerated from (3 = 0.956 (7 - 3.4) to j3 = 0.995 (7 - 10.3). 
The increase of velocity is less at smaller viewing angles. 

An alternative model for the acceleration and brightening of 
component E would be a change of the jet inclination to the 
line of sight from about 24° to about 1 1° at this l ocation in the 
jet as observed in the case of the quasar 3C 279 (Hom an et alj 
Bool . However, we see no significant change in jet position 
angle in the sky which would be expected to accompany such 
a large change in jet direction. Moreover, subsequent com- 
ponents, particularly G, do not show the same kind of large 
acceleration in this region. 

Direct identification of component E with component El 
is not straightforward in the frame of expansion, as compo- 
nent El in epoch 1999.38 was smaller than component E in 
1998.18 (see Table 2). However, component E3 in epoch 
1999.38 is larger than component E in 1998.18. We can in- 
terpret this as component E including components El and E3 
(and maybe E4). These components would be i ndistinguish- 
able in our observations before 1999.38. In fact, Jorst ad et alj 
(2005) monitored 3C111 between 1998 and 2001 with the 
VLBA at 43GHz. They find an emission complex, that can 
be identified with our component E, that gradually stretches 
out as it travels from ~ 2 mas from the core in 1998 to 
roughly between 5 mas and 8 mas from the core in 2001. 
Their leading component CI can be identified with our com- 
ponent El, their component c2 with E2 and their cl with 
E3. At their higher angular resolution, Jorstad et al. can 
separate components C 1 and c 1 already in early 1998. In 
agreement with our analysis at 15 GHz, they detect c2 (E2) 
about a year after they detect cl (E3). They do not detect 
a component corresponding to E4 but this may be an effect 
of partially resolving out the jet structure at their higher ob- 
serving frequency, particularly in later epochs. It is further 
interesting to note that the observed speeds at both frequen- 
cies agree well. ForEl(Cl), /i a pp,2cm = 1.69 ± 0.04 mas yr" 1 
at 15 GHz and jU a pp,7min = 1. 77 ±0.06 mas yr" 1 at 43 GHz; 
for E2(c2), /i a pp,2cm = 1.29 ± 0.06 mas yr" 1 at 15 GHz and 
Aiap P ,7mm = 1 -23 ± 0.04 mas yr" 1 at 43 GHz; and for E3(cl), 
Map P ,2cm = 1.22 ± 0.05 mas yr" 1 at 15 GHz and /i ap p,7mm = 
1. 07 ±0.02 mas yr" 1 at 43 GHz. The discrepancy in the speeds 
measured for E3 and c 1 seems to be due to a slight accelera- 
tion of E3 after 2002. A fit to the 15 GHz data of E3 between 
1999 and 2002 alone yields a slower speed of ~ 1.0 mas yr" 1 
similar to the speed of cl in the same time period at 43 GHz. 
In their work, Jorstad et al. do not report acceleration of com- 
ponents from 2 mas to 4 mas. However, this is likely due to the 
fact that their observations started in early 1998, thus missing 
the first observations of component E presented in this paper, 
when its speed has been measured to be smaller. 

4.3. Between 2 pc and 6 pc: Recollimation of the Jet 

Inspection of Fig. 9 shows that the back-extrapolation of 
the brightness temperature of component E from regime II to 
regime I is at least two orders of magnitude too high if this 
extrapolation is based on the gradient given by El. The low 
brightness temperature of component E in regime I cannot be 
explained by opacity effects because the radio-light curve in 
Fig. 1 shows that the source was optically thin from 1997 on. 
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Moreover, if we identify component E with E 1 , it is Doppler- 
deboosted from epoch 1998. 18 to 1999.38 due to the accelera- 
tion and a relatively large viewing angle; thus, we are not able 
to explain the increase in brightness temperature in terms of 
Doppler boosting. However, compact sub-components may 
have larger brightness temperatures, so that the 7b values plot- 
ted in Fig. 9 for E in regime I inward of about 3 mas may 
represent lower limits for compact components already em- 
bedded in the unresolved structure. 

Not only E/El but also components G and H show 
an increase in total flux density several milliarcseconds 
downstream. Compared to E/El, these somewhat weaker 
components exhibit their flux-density maxima at somewhat 
larger distances from the core (compare Fig. 6 and Fig. 7). 
This can be explained if the gas in the components trav- 
els through a mild standing shock in a recollimation re- 
gion. This eff ect has been observ ed in numerical simula- 
tions of parsec (|G6mez et a"Tl ll997) and kiloparsec scale jets 
dPerucho & MartHl2007l) . The material in the components is 
expected to be overpressured with respect to its environment, 
thus expanding into it. After the initial expansion, the com- 
ponents become underpressured with respect to the underly- 
ing flow. The resulting recollimation leads to the formation 
of a shock, whose strength depends on the initial degree of 
overpressure of the material in the component. This process 
explains the increase in flux density and brightness tempera- 
ture as due to compression of the gas in the recollimation. In 
Figs. 6 and 7, we see that the flux density of component E in- 
creases closer to the core than for component G and H, which 
is consistent with the for mer being slower than th e latter, thus 
recollimating earlier (see Per ucho & Martil l2007). It also ex- 
plains why we see a significant acceleration only in the faster 
expanding, brighter component E/El . 

Finally, after this mild recollimation, the fluid becomes 
overpressured with respect to its environment, thus further ex- 
panding and accelerating downstream. 

4.4. Near 3 pc: The Role of the External Medium 

The polarization behavior of components H and I can be 
understood in terms of an interaction between the jet and the 
external medium at a distance of 3.3 mas (= 3.3 pc) in the jet. 
Assuming no Faraday rotation, the EVPA of component H 
within approximately 3.3 mas from the core indicates a trans- 
verse magnetic field order as might be expected for a trans- 
verse shock propagating down the jet. The change in the 
fractional polarization, its north-south gradient, and the ro- 
tation of the EVPA suggest that a contact surface persists at 
the southern boundary of the jet beam at a distance of ap- 
proximately 3.3 mas downstream the jet core. If the bulk jet 
material flows faster than the flow at the southern boundary, 
the magnetic field is stretched through shear. Our overall pic- 
ture then is of an originally transverse shock interacting with 
the jet on the southern side of the jet at 3.3 mas from the 
core. The interaction changes the component's magnetic field 
through some combination of oblique shock and differential 
flow resulting in a magnetic field approximately parallel to 
the jet axis in the later epochs. No strong shock is needed at 
this location in the jet but this reg ion m ay be identified with 
the recollimation region (see Sect. 14.3b at about the same po- 
sition in the jet). In this picture, the jet-medium interaction 
may form an effective nozzle which accelerates the jet on one 
edge relative to the other. 

An alternative explanation for the observed polarization 
structure and dynamics of 3C 111 can be found by consid- 



ering an inhomogeneous external Faraday screen. Such a 
screen could produce the observed differential rotation of 
the EVPA wh ile a component travels through a given region 
along the jet. IZavala & Tavlori d2002l) observed 3C 1 1 1 with 
the VLBA and produced a Faraday rotation-measure map be- 
tween 8 GHz and 15 GHz. They find strong Faraday rotation, 
~ 730radm~ 2 , at the same distance from the core (3.3 mas) 
where our observations show the swing of the EVPA of the 
component H and steeply decreasing Faraday rotation further 
downstream. However, we note that 730radm" 2 translates to 
17° of rotation at 15 GHz which alone is not enough to ex- 
plain the change in EVPA that we observe while component 
H travels through this region. On the other hand, the steep de- 
crease of the Far aday rotation measured up- and downstream 
of this region bv IZavala & Tavlori d2002l) again agrees with a 
change of the external gas density at this point, which in turn 
may be identified with the pressure gradient responsible for 
the component expansion and acceleration. 

A combination of inhomogeneous Faraday rotation and an 
interaction between the jet plasma and its ambient medium 
appears most likely to explain our observations of the vary- 
ing linear polarization structure; however, both explanations 
point to the role of the external medium, either through a dis- 
crete interaction or a rapid decrease in external gas pressure, 
in shaping the jet flow downstream of this location. 

4.5. Between 3 pc and 5 pc: Formation of Trailing 
Components 

The components E2, E3, and E4 can be interpreted as trail- 
ing components forming in the wake of the leading El which 
is identical with the original component E. This scenario is 
attractive becaus e the basic concep t of trailing components 
as introduced by lAgudo et al.l d2001l) predicts the formation 
of trailing features in the wake of the initial perturbation in 
the jet flow. Such a behaviour has first been found both as- 
sociated with bright sub- a nd superluminal jet components in 
Centaurus A and 3C 120 (jTingay. Preston. & Jaunceyl 120011 : 
Gomez et al. 200l])- iJorstad et al.1 (j2005) report trailing com- 
ponents in four additional sources (3C273, 3C345, CTA 102 
and 3C 454.3) and in 3C 111 (see below). 

The interaction of the external medium with a strong shock 
pinches the surface of the jet, leading to the production of 
pinch-body mode Kelvin- Helmholtz instabilities: the trailing 
features. Hence, a single strong superluminal shock ejection 
from the jet nozzle may lead to the production of a multiple 
set of emission features through this mechanism. The trailing 
features have a characteristic set of properties, which make 
them recognizable with high resolution VLBI: they form in 
the wake of strong components instead of being ejected from 
the core of VLBI jets, they are related to oblique shocks, they 
are always slower than the leading feature, and (if the under- 
lying jet has a certain opening angle) they should be generated 
with a wide range of apparent speeds (from almost stationary 
near the core to sup erluminal further downstream). Moreover, 
Agud cTeTail d2001l) showed that the separation between the 
trailing components increases downstream due to their mo- 
tion down a pressure gradient. 

All this is in agreement with what we observe in the trail- 
ing components of E 1 and with our interpretation of an ex- 
pansion of the jet in a density decreasing ambient medium. 
For the tim e range covered by their observations (1998.23 to 
2001 .28t. ljorstad et al.1 d2005l) also identified the trailing phe- 
nomenology in this source. 

The north-south gradients detected in the linear- 
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polarization emission in the region where the trailing 
features are formed, is in agreement with an oblique shock 
structure. The steep brightness-temperature gradients of the 
trailing components indicate that the particle and magnetic 
field density associated with these components evolve in 
a different way compared to the "normal" jet flow. These 
shocked regions may be more overpressured with respect 
to their environment, making them expand rapidly. This 
fast expansion implies a larger positive value of /, which, 
however, is compensated by an even larger (negative) value 
of n and b(l-a) in equation Q] resulting in a very steep 
brightness temperature gradient (regime II). 

Pinching modes of the Kelvin-Helmholtz instability were 
shown to coupl e to the trail i ng co mponents observed in the 
simulations in Agu do et al.l d2001). In the case of compo- 
nents E2-E4, the distance between them ranges from 0.7-0.8 
mas at the first epochs in which they are observed, to almost 
2.0 mas in the latest epochs. Taking into account that: a) 
their FWHM is of the same order (see Table 2); b) that these 
wavelengths have to be corrected for geometrical and rela- 
tivistic effects, resulting in a maximum intrinsic wavelength 
of ~ 0.7 mas, and c) that the size of the c omponents can be of 
the or der or smaller than the jet radius (Perucho & Lobanov 
120071) . this implies coupling of the pinch ing to wavelengths of 
the order or smaller than the jet radius. Perucho et al. (2007) 
have shown that resonant Kelvin-Helmholtz instabilities as- 
sociated to high-order body modes appear in sheared jets at 
these wavelengths. These modes have larger growth rates than 
low-order body modes or surface modes, and their growth 
brings the jet to a final quasi-steady state in which it remains 
well-collimated and generates a hot shear-layer which shields 
the core of the jet from the ambient medium. Interestingly, the 
jet in 3C 1 1 1 is known to be well-collimated up to kiloparsec 
scales. Further research in this direction is needed in order 
to check the influence of the resonant modes in the long term 
evolution of this jet. 

A by-product of the interpretation of these components as 
Kelvin-Helmholtz instabilities is the fact that it allows us to 
put constraints to the velocity of the jet. We can regard 
the wave speed as the minimum speed of the jet flow, as 
KH modes have an upper limit in their wave speeds that is 
precisely the veloci ty of the flow in which they propagate 
dPerucho et al.l20 06). The upper limit is given by the speed of 
El, interpreted as a shock wave, that has to be thus faster than 
the underlying flow. In this picture, we would have the struc- 
ture El moving with Lorentz factor 7 ^ 8.3 through a jet with 
Lorentz factor 8.3 > jj > 4.6 in the accelerated region (post 
1999.38), where the lower limit is given by the Lorentz factor 
of component E2, the fastest of the three trailing components 
identified here. 

5. SUMMARY AND CONCLUSIONS 

In this paper, we have investigated the parsec-scale jet kine- 
matics and the interaction of the jet with its ambient medium 
in the broad-line radio galaxy 3C 111. Our analysis has 
demonstrated that a variety of processes influence the jet dy- 
namics in this source: a plasma injection into the jet beam 
associated with a major flux-density outburst leads to the for- 
mation of multiple shocks that travel at different speeds down- 
stream and interact with each other and with the ambient 
medium. The primary perturbation causes the formation of 
a forward and a backward shock (or rarefaction). The latter 
fades away so fast that is likely to remain undetected in minor 
ejections. A separate work by iPerucho et al.l d2007l) focuses 



on the nature and characteristics of these initial components. 
Several parsecs downstream, the jet plasma enters a region of 
rapidly decreasing external pressure, expands into the jet am- 
bient medium and accelerates. In the following, the plasma 
gets recollimated and trailing features are formed in the wake 
of the leading component. 

A particularly interesting aspect of the source 3C 111 
in the light of this and other recent works is that it is 
one of the very rare non-blazar gamma-ray bright AGN. 
Besides Centaurus A dSreekumar et all 1 19991) and the pos- 
sible identificatio n of NGC6251 with the EGRET source 
3EG J 162 1+8203 dMukherjee et alJl2002h . 3C 111 is the only 
AGN whose jet-system is inclined at a relatively large angle 
to the line of sight and th at has a reliable EGRET identifica- 
tion: ISguera et al. (2005) reconsidered the possible identifi- 
cation of the EGRET source 3EGJ04 16+3650 wit h 3C 111, 
which was first suggested by Hartman et al.l (1 19991) but con- 
sidered unlikely because of the poor positional coincidence. 
Very recently, R. C. Hartman & M. Kadler (in prep.) found 
that 3EG J0416+3650 is composed out of at least two dis- 
tinct components. One of them is the dominant source above 
1 GeV and is in excellent positional agreement with the loca- 
tion of 3C 1 1 1 . Compared to blazars, the large inclination an- 
gle and the relatively small distance of 3C 1 1 1 allow us to re- 
solve structures along the jet that are as small as parsecs in de- 
projection and which would be heavily blended with adjacent 
features in blazar jets. As demonstrated in this paper, VLBA 
observations of 3C 111 probe a variety of physically differ- 
ent regions in a relativistic extragalactic jet such as a com- 
pact core, superluminal jet components, recollimation shocks 
and regions of interaction between the jet and its surrounding 
medium, which are all possible sites of gamma-ray produc- 
tion. From early 2008 on, the gamma-ray satellite GLAST 
(lLott et al.l [2007) is going to monitor the sky. If detected by 
GLAST, 3C 111 may become a key source in the quest for an 
understanding of the origin of gamma-rays from extragalac- 
tic jets. In addition, the combination of GLAST and VLBA 
data with spectral data at intermediate wavelengths (optical, 
IR, X-ray) may allow a better determination of jet parameters 
and relativistic beaming effects than in most blazars because 
of the higher linear resolution offered by this nearby and only 
weakly projected jet system. 

Our observations of 3C 1 1 1 are qualitatively in remarkable 
agreement with numerical relativistic hydrodynamic struc- 
tural and e mission sim u lation s of j ets such as the on es pre- 
sented by lAgudo et all d2001l) and lAlov et al.l d2003l) . Fur- 
ther progress is being made in the transition from two- 
dimensional to three-dimensional simulations of relativistic 
jets and in the development of new methods considering 
magnetic fields (fLeismann et alJ 120051 : iMizuno et alJ 120071 : 
iRoca-Sogorb et aHl2008l e.g.,), the e quation of state for rela- 
tivistic gases (Perucho & Marti 2007), and radiative processes 
(e.g.. iMimica et al. 2004, 2007, and Mimica et al. in prepa- 
ration). But so far neither observational data nor simulations 
have reached an adequate level of detail and completeness in 
order to allow us a quantitative direct comparison of numer- 
ical models and observed relativistic jet structure and evolu- 
tion. In particular, it is not feasible today to fit iteratively the 
parameters of relativistic magneto-hydro-dynamical (RMHD) 
jet simulations to match the brightness distribution observed 
for any individual source. The main reasons for this are a) 
the immense computational power required to conduct a re- 
alistic (i.e., sufficiently detailed) modern 3D jet simulation 
and b) the highly non-linear nature of RMHD plasmas and 
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their evolution. Simulation results depend critically on the 
starting conditions like the exact velocity, composition, and 
transversal structure of the flow, the structure and strength of 
the magnetic field and the jet environment. Future develop- 
ment of computational power will allow us to use larger res- 
olutions to decrease the numerical viscosities, and to imple- 
ment nonlinear and microphysics processes into simulations. 
VLBA observations are capable of putting hard quantitative 
constraints on the input parameters for RMHD jet simulations 
if they are densely sampled over several years. Polarimet- 
ric observations at multiple radio frequencies may allow the 
effects of jet-intrinsic magnetic-field variations and external 
Faraday-screen inhomogenities or temporal variations to be 
disentangled. Such data at 15 GHz are on the way, e.g., as part 
of the next phase of the MOJAVE program, in which rapidly 
evolving sources like 3C 111 are being observed every two 
months. 
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APPENDIX 

A. THE JET INCLINATION ANGLE 

The 1996 radio outburst of 3C 1 1 1 puts strong constraints on the angle to the line of sight for this source, if one assumes that a 
similarly bright component as E has been ejected in the counterjet, as well. Due to differential Doppler boosting, the flux density 
ratio between the jet- and counter-jet emission is 



S } _ /l+(3cos9 
Scj \l-/3cos0 



Thus, for a given jet to counter-jet ratio x = 



/3cos(9: 



x-l 
x+ 1 



(Al) 



(A2) 



With a = 0.3, Sj = 3.4 Jy (components E and F in 1997.19), Scj < lOmJy and /3 < 1, 9 < 21°. For a realistic jet speed of, e.g., 
(3 = 0.956 (7 = 3.4), the angl e to the line of sight is: 9 = 19°. An estimate close to this value can be derived from the variability 
Doppler factor me asured by [Lahteenmaki & Valtaoial ( 1 1999b and the apparent superluminal jet speed. As outlined in detail in 
ICohen et al.l (120071) . this leads to a value of 9 ~ 15°. 

It is important to note that this calculation implicitly assumes symmetry between the jet and counter-jet, which in projection 
does not have to be the case if the cou nter-jet is covered by an obscuring torus as it is well-established for systems at larger 
incl ination angles (e.g., NG C 1052: see Ka dler et al.l l2004). Indeed, Faraday rotati on measuremen ts towards the 3C 111 pc-scale 
jet dZavala & Tavlorir2002t see also Sect. |4.4| > and X-ray spectral observations dLewis et al.l 120051) suggest substantial amounts 
of obscuring material. Free-free absorption could also substantially lower the counter-jet radio emission and allow for larger jet 
angles to the line of sight. 

An independent lower limit on the inclination angle of 9 > 21° was given by lLewis et all (120051) assuming that the deprojected 
size of the largescale 3C 111 double-lobe structure is smaller than 500/T 1 kpc. This discrepancy implies that either 3C 111 is 
unusually large or there is a misalignement between the large-scale jet-axis and the parsec-scale jet axis inclination to the line of 
sight, although the projected position angles of the large-scale jet (63°) and the parsec-scale jet (~ 65°) are almost the same. 
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B. IMAGE-PLANE VS. (t/,V)-PLANE MODEL FITTING 

It is interesting to compar e our results from this v ery detailed analysis of one individual object with the results of the 
kinematical-survey analysis of Kellerm ann et al.l d2004l) . who investi gated the speeds of 1 10 extragalactic jets, including 3C 111, 
based on the VLBA 2 cm Survey data between 1994 and 2001. Kellerman n et alJ (120041) made the component identification 
in the image plane and represented the evolving structure of component E and its trailing features by only one component 
whereas, in this work w e distinguish the sub-components E, F, El, E2, E3, and E4. Formally, the speed of (4.9 ± 0.2)c found by 
iKellermann et alJ d2004h is in good agreement with the speed of our leading component, El, so the much simpler model derived 
from image-plane analysis represents the fastest moving structure. The acceleration (with respect to the ejecta's smaller velocity 
prior to mid 1999), as well as the additional components that we interpret as a backward shock and trailing jet features become 
visible only after a more complicated model fitting of the data in the (u, v)-domain. The necessarily less-complex model used in 
a survey analysis like the one conducted by Kellerma nn et alJ ((2004) is only part of the reason for this discrepancy. Image-plane 
analysis makes it very difficult to interpret a moving feature that changes its structure in a complex way and that has no clear 
persistent brightness maxima. In addition, {u, v)-plane fitting in general achieves higher angular resolution so that the two bright, 
but closely separated components E and F could not be distinguished in in early to mid 1996. Figure [TT| demonstrates that even 
in October 1996 when E and F are separated by only 0.3 mas and are located within 1 mas from the core, a one-component model 
clearly fails to represent this compact structure. 



Visibility Amplitude 




100 200 300 400 

Radial UV Distance Along P.A. 70 deg 
[Million Wavelengths] 



100 200 300 400 

Radial UV Distance Along P.A. 70 deg 
[Million Wavelengths] 



FIG. 1 1 . — Visibility amplitudes projected radially along P.A. 70° for the 1996.82 observation of 3C 111. The double-peak indicates the presence of a bright 
structure within less than a milliarcsecond of the core. The top left panel shows a model (black) which was fitted to the data (green) consisting of one model 
component for the newly ejected jet feature and the residuals of this model are shown in the bottom left panel. Up to about 600 mjy of correlated flux density is 
missing from the model. The right panels show the same data fitted by a model consisting of two components (corresponding to E and F). 
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TABLE 1 

Journal of VLBA 2cm Survey observations of 3C 1 1 1 analyzed in this study. 



Epoch 


Code 


5tot 


rms 


C* 


m* 






[Jy] 


[mJy/beam] 


[mJy/beam] 


[%] 


1995.27 


BK016 


2.6 


1.8 


6.5 




1995.96 3 


BK037A 


1.8 


1.1 


3.4 




1996.82 


BK37D 


4.8 


0.5 


2.1 




1997.19 


BK048 


6.0 


0.9 


4.9 




1997.66 


BK052A 


4.1 


1.6 


6.4 




1998.18 


BK052B 


2.0 


2.1 


5.5 




1999.39 


BK068A 


2.6 


0.6 


1.9 




1999.85 


BK068C 


2.4 


1.2 


4.1 




2000.49 


BT051 


2.6 


0.5 


1.4 


1.5 


2001.17 


BK068E 


2.1 


0.3 


1.3 




2002.00 


BR077D 


2.2 


0.3 


1.1 




2002.19 


BR 0771 


2.2 


0.2 


0.9 




2002.77 


BL111C 


1.6 


0.5 


1.7 


0.6 


2003.65 


BL111J 


2.1 


0.3 


1.2 


0.3 


2004.27 a 


BL111L 


1.9 


0.4 


1.4 


0.7 


2004. 80 b 


BL111P 


2.7 


0.3 


1.0 


0.7 


2005.02 


BL 123A 


3.0 


0.5 


1.8 


0.5 


2005.73 


BL 1230 


2.4 


0.3 


1.2 


0.7 



* Lowest contour in Fig. [3] 
t Degree of polarization (see Fig. [Tol 
a No data from antenna at Mauna Kea 
b No data from antenna at St.Croix 
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TABLE 2 
Model fit parameters 



ID 


Flux Density [mJy] 


Radius [mas] 


P.A. a [°] 


FWHM [mas] 


ratio 


4>n 


1995 77 





1371.04 ± 205.66 


0.00 


24.00 


0.33 ± 0.60 


0.43 


56.22 


1 1 
1 > 


n£ 4- 111 A 1 


0.61 ± 0.30 


61.18 


mi _i_ a c/i 

U.Zl It U.J4 


1 nn 




c 


149 SQ + ?2 44 


1.14 ±0.21 


60.80 


41 + () 64 


1 00 




B 


112.85 ± 16.93 


2.06 ± 0.29 


69.41 


0.62 ± 0.79 


too 




A 


73.30 ± 10.99 


3.97 ± 0.30 


70.69 


0.00 ± 0.50 


1.00 


— 






1995 96 











1173.64 ± 176.05 


0.00 




0.35 ± 0.12 


0.59 


47.42 


D 


353.66 ± 53.05 


0.73 ± 0.30 


61.85 


0.30 ± 0.12 


l b 




c 


145.60 ± 21.84 


1.72 ±0.21 


62.73 


0.38 ± 0.13 


l b 




B 


46.15 ± 6.92 


3.37 ± 0.29 


70.25 


0.61 ± 0.16 


I b 




A 


43.54 ± 6.53 


4.53 ± 0.30 


69.99 


_c 


l b 








1996 87 











1645 03 ± 246 75 


0.00 




26 ± 1 1 


17 


60 63 


F 


1608.94 ± 241.34 


0.44 ± 0.07 


54.71 


0.31 ± 0.12 


0.19 


87.46 


E 


1344.91 ± 201.74 


0.73 ± 0.02 


67.83 


0.23 ±0.11 


0.61 


35.24 


XI 


121.29 ± 18.19 


1.06 ±0.3 


61.59 


0.33 ± 0.12 


l b 


- 


C 


O O l*\ A 1 CIA 

38.24 ± 5.74 


3.28 ± 0.21 


66.12 


0.88 ± 0.20 


ib 


— 


B 


38.04 ± 5.71 


4.76 ± 0.06 


68.08 


0.66 ± 0.16 


l b 


— 






I9Q7 19 











2490.20 ± 373.53 


0.00 




0.38 ±0.13 


0.26 


62.36 


F 


1980.30 ± 297.05 


0.77 ± 0.07 


55.51 


0.34 ±0.12 


0.39 


57.52 


E 


1410.05 ± 211.51 


1.06 ± 0.02 


71.09 


0.34 ± 0.12 


0.35 


51.00 


C 


32.15 ±4.82 


3.72 ± 0.21 


65.48 


0.23 ±0.11 


l b 




B 


35.26 ± 5.29 


5.45 ± 0.06 


67.95 


0.72 ±0.18 


l b 








1997 66 











1697.09 ± 254.56 


0.00 




0.51 ± 0.14 


C 


59.71 


F 


940.20 ± 141.03 


0.93 ± 0.07 


60.22 


0.27 ±0.11 


l b 




E 


1461.49 ± 219.22 


1.56 ± 0.02 


67.18 


0.40 ± 0.13 


l b 








1998 18 











1128.63 ± 169.29 


0.00 




0.82 ±0.19 


0.47 


6.27 


F 


131.71 ± 19.76 


1.36 ±0.07 


59.83 


0.59 ± 0.16 


l b 




E 


815.70 ± 122.35 


2.07 ± 0.02 


67.53 


0.90 ± 0.21 


0.65 


1L93 






1999 19 












1509.63 ± 226.44 


0.00 




0.38 ± 0.13 


C 


64.00 


X3 


225.98 ± 33.90 


0.59 ± 0.30 


62.46 


0.14 ±0.10 


l b 




X2 


54.38 ± 8.16 


1.25 ± 0.30 


65.08 


0.56 i 


i 0.15 


l b 




E3 


470.42 ± 70.56 


3.51 ± 0.28 


68.80 


1.34 i 


i 0.29 


0.37 


57741 


El 


298.52 ± 44.78 


4.21 ± 0.12 


66.02 


0.56 i 


i 0.15 


0.42 


35.72 



.1999.85. 






960.22 ± 144.03 


0.00 




0.78 ± 0.19 


0.41 


-4.69 


X4 


743.68 ±111.55 


0.55 ± 0.30 


59.07 


1.13 i 


i 0.25 


0.58 


61.59 


E3 


349.03 i 


i 52.35 


3.82 i 


i 0.28 


67.78 


1.26 i 


i 0.27 


0.69 


50.04 


E2 


122.27 i 


i 18.34 


4.66 i 


i 0.21 


65.12 


0.84 i 


i 0.20 


0.39 


-16.32 


El 


267.68 i 


i 40.15 


5.12 i 


i 0.12 


64.17 


0.84 - 


i 0.19 


0.40 


-7.05 










?nnn 49 













1668.98 ± 250.35 


0.00 




0.33 ± 0.12 


0.11 


66.21 


X7 


249.65 ± 37.45 


0.65 ± 0.30 


56.02 


0.45 ±0.13 


l b 




X6 


133.82 ± 20.07 


1.31 ± 0.30 


59.94 


0.48 ±0.14 


l b 




X5 


43.07 ± 6.46 


2.34 ± 0.30 


66.11 


0.71 ± 0.17 


l b 




E4 


66.54 ± 9.98 


3.80 ± 0.15 


68.04 


0.47 ±0.14 


l b 




E3 


160.55 ± 24.08 


4.46 ± 0.28 


63.86 


0.61 ± 0.16 


l b 




E2 


125.58 ± 18.84 


5.35 ± 0.21 


64.62 


0.74 ± 0.18 


l b 




El 


135.77 ± 20.37 


6.08 ±0.12 


64.13 


0.68 ±0.17 


l b 





^001.17. 






1596.68 ± 239.50 


0.00 




0.44 ±0.13 


0° 


61.41 


G 


243.56 ± 36.53 


1.13 ± 0.07 


58.96 


0.88 ± 0.20 


0.22 


54.85 


E4 


91.23 ± 13.68 


3.96 ± 0.15 


66.28 


1.11 ± 0.24 


0.35 


48.34 


E3 


99.76 ± 14.96 


5.21 ± 0.28 


62.80 


1.24 ± 0.27 


0.54 


65.22 


E2 


45.34 ± 6.80 


6.20 ± 0.21 


66.06 


0.91 ± 0.21 


0.92 


55.42 


El 


45.18 ± 6.78 


7.20 ± 0.12 


64.16 


0.86 ± 0.20 


0.83 


75.97 



J2002.00. 






1514.03 ± 227.10 


0.00 




0.49 ±0.14 


0.09 


61.18 


H 


424.80 ± 63.72 


1.02 ± 0.02 


64.67 


0.46 ± 0.14 


0.39 


66.86 


G 


60.57 ± 9.09 


2.59 ± 0.07 


58.88 


0.94 ± 0.21 


0.33 


50.61 


E4 


121.10 ± 18.16 


4.72 ± 0.15 


64.66 


1.40 ± 0.30 


0.35 


47.76 


E3 


19.67 ± 2.95 


5.90 ± 0.28 


65.92 


0.82 ±0.19 


l b 




E2 


21.62 ±3.24 


6.98 ± 0.21 


64.03 


1.12 ±0.25 


l b 




El 


24.46 ± 3.67 


8.43 ±0.12 


65.14 


1.04 ± 0.23 


l b 





J.002.19. 






1688.54 ± 253.28 


0.00 




0.56 ±0.15 


0.09 


63.02 


H 


317.77 ±47.66 


1.31 ± 0.02 


64.52 


0.48 ± 0.14 


0.56 


65.20 


G 


29.58 ± 4.44 


3.07 ± 0.07 


57.88 


_C 


l b 




E4 


118.14 ± 17.72 


4.91 ± 0.15 


64.22 


1.51 ±0.32 


0.36 


52?47 


E3 


12.28 ± 1.84 


6.18 ± 0.28 


67.59 


0.51 ± 0.14 


l b 




E2 


17.46 ± 2.62 


7.77 ± 0.21 


64.31 


1.00 ±0.22 


l b 




El 


12.35 ± 1.85 


8.74 ±0.12 


66.18 


0.54 ±0.15 


l b 





J.002.77 . 






1280.01 ± 192.00 


0.00 




0.47 ±0.14 


C 


62.23 


J 


100.40 ± 15.06 


0.61 ± 0.28 


68.62 


0.30 ± 0.12 


l b 




I 


78.58 ± 11.79 


1.30 ± 0.22 


61.52 


0.34 ± 0.12 


l b 




H 


90.33 ± 13.55 


2.30 ± 0.02 


64.39 


0.50 ± 0.14 


l b 




G 


47.49 ± 7.12 


4.16 ±0.07 


61.40 


0.84 ± 0.20 


l b 




E4 


53.64 ± 8.05 


5.51 ± 0.15 


63.49 


0.87 ± 0.20 


l b 




E3 


13.63 ± 2.04 


6.85 ± 0.28 


64.21 


0.84 ± 0.20 


l b 
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ID Flux Density [mJy] Radius [mas] P.A. a [°] FWHM [mas] ratio </>[°] 



E2 


1 84 + 1 63 


8.51 ±0.21 


66.27 


1 23 ± 27 


l b 




El 


11314-1 70 

A L .J A _1_ 1 . / \J 


9.69 ±0.12 


62.75 


1 45 4- n 3 1 

J. _i_ u.-ji 


l b 



















1046.67 ± 157.00 


0.00 




0.32 ±0.12 


0.20 


74.73 


X9 


576.24 ± 86.44 


0.36 ± 0.30 


68.30 


0.09 ±0.10 


l b 




K 


177.19 ± 26.58 


0.69 ± 0.20 


66.00 


0.14 ± 0.10 


l b 




X8 


66.75 ± 10.01 


1.30 ± 0.30 


62.77 


0.17 ±0.11 


l b 




J 


34.80 ± 5.22 


1.94 ± 0.28 


63.40 


0.33 ± 0.12 


l b 




I 


21.58 ±3.24 


2.71 ± 0.22 


62.07 


0.51 ±0.14 


l b 




H 


67.76 ± 10.16 


3.86 ± 0.02 


62.83 


0.97 ± 0.22 


l b 




G 


64.55 ± 9.68 


5.73 ± 0.07 


60.73 


0.59 ±0.15 


l b 




E4 


7.48 ± 1.12 


6.46 ±0.15 


60.88 


0.57 ±0.15 


l b 




E3 


13.78 ± 2.07 


7.95 ± 0.28 


65.78 


1.59 ± 0.33 


l b 




El 


11.72 ± 1.76 


10.90 ±0.12 


63.95 


1.79 ±0.37 


l b 





2004.27 






1204.12 ± 180.62 


0.00 




0.31 ±0.12 


0.14 


70.05 


X10 


270.09 ± 40.51 


0.44 ± 0.30 


68.51 


_C 


l b 




L 


105.50 ± 15.83 


0.82 ± 0.20 


69.58 


0.20 ± 0.11 


l b 




K 


38.46 ± 5.77 


1.58 ± 0.20 


66.52 


0.41 ±0.13 


l b 




J 


26.88 ± 4.03 


2.42 ± 0.28 


63.12 


0.40 ±0.13 


l b 




I 


11.76 ± 1.76 


3.36 ± 0.22 


65.95 


1.09 ± 0.24 


l b 




H3 


110.03 ± 16.50 


4.45 ± 0.21 


63.91 


0.45 ± 0.14 


l b 




H2 


76.83 ± 11.53 


4.74 ± 0.18 


59.02 


0.35 ±0.12 


l b 




HI 


26.80 ± 4.02 


5.45 ± 0.24 


62.26 


0.82 ±0.19 


l b 




G 


26.20 ± 3.93 


6.81 ± 0.07 


61.25 


1.12 ± 0.25 


l b 




E4 


15.80 ±2.37 


6.79 ±0.15 


61.23 


0.58 ±0.15 


l b 




E3 


5.15 ±0.77 


8.63 ± 0.28 


63.85 


0.75 ±0.18 


l b 




E2 


5.23 ± 0.78 


10.21 ± 0.21 


67.04 


1.37 ±0.29 


l b 




El 


3.92 ± 0.59 


12.94 ± 0.46 


64.45 


0.77 ±0.18 


l b 





2004.80 






1590.66 ± 238.60 


0.00 




0.28 ±0.11 


0.38 


66.10 


N 


507.49 ±76.12 


0.33 ± 0.15 


61.95 


0.14 ± 0.10 


l b 




M 


294.28 ±44.14 


0.56 ± 0.15 


70.23 


0.22 ±0.11 


l b 




L 


34.00 ±5.10 


1.06 ± 0.20 


65.05 


0.30 ± 0.12 


l b 




K 


15.62 ±2.34 


1.75 ± 0.20 


70.84 


0.38 ± 0.13 


l b 




J 


13.76 ± 2.06 


2.55 ± 0.28 


66.12 


0.33 ±0.12 


l b 




I 


20.41 ± 3.06 


3.62 ± 0.22 


64.10 


0.62 ±0.16 


l b 




H4 


41.27 ± 6.19 


4.93 ± 0.24 


63.72 


0.64 ± 0.16 


l b 




H3 


88.46 ± 13.27 


5.25 ± 0.21 


62.60 


0.87 ± 0.20 


l b 




H2 


97.12 ± 14.57 


5.59 ±0.18 


59.49 


1.10 ± 0.24 


l b 




HI 


13.27 ± 1.99 


6.47 ± 0.24 


61.41 


1.82 ±0.38 


l b 




G 


21.67 ± 3.25 


7.80 ± 0.07 


61.58 


0.37 ±0.12 


l b 




E3 


4.32 ± 0.65 


9.85 ± 0.28 


66.41 


0.97 ± 0.22 


l b 




El 


6.74 ± 1.01 


13.19 ±0.46 


64.15 


0.77 ±0.18 


l b 





2005.02 






1631.40 ± 244.71 


0.00 




0.24 ±0.11 


0.21 


61.51 


N 


542.58 ± 81.39 


0.31 ± 0.15 


70.03 


0.06 ±0.10 


l b 




M 


386.90 ± 58.03 


0.67 ±0.15 


63.61 


0.15 ± 0.10 


l b 




Xll 


187.02 ± 28.05 


0.86 ± 0.30 


70.85 


0.20 ±0.11 


l b 




L 


17.02 ± 2.55 


1.47 ±0.20 


67.37 




l b 




K 


12.09 ± 1.81 


1.99 ±0.20 


68.81 


0.39 ± 0.13 


l b 




J 


16.23 ± 2.43 


3.21 ± 0.28 


64.59 


0.99 ± 0.22 


l b 




I 


12.46 ± 1.87 


4.28 ± 0.22 


64.94 


0.54 ±0.15 


l b 




H4 


101.42 ± 15.21 


5.42 ± 0.24 


62.82 


0.47 ±0.14 


l b 




H3 


58.52 ± 8.78 


5.81 ± 0.21 


59.93 


0.39 ±0.13 


l b 




H2 


46.72 ± 7.01 


6.18 ±0.18 


59.53 


0.79 ± 0.19 


l b 




HI 


4.59 ± 0.69 


7.16 ±0.24 


61.98 


_C 


l b 




G 


15.88 ± 2.38 


8.07 ± 0.07 


61.54 


0.87 ± 0.20 


l b 




E3 


4.09 ± 0.61 


10.02 ± 0.28 


66.09 


1.51 ± 0.32 


l b 




El 


5.04 ± 0.76 


13.35 ± 0.46 


65.18 


1.28 ±0.27 


l b 








7005 7? 











1197.32 ± 179.60 


0.00 




0.28 ±0.11 


0.30 


69.16 


X14 


480.43 ± 72.06 


0.35 ± 0.30 


70.12 


0.12 ±0.10 


l b 




X13 


224.64 ± 33.70 


0.67 ± 0.30 


68.78 


0.13 ±0.10 


l b 




N 


71.31 ± 10.70 


1.07 ± 0.15 


63.01 


0.23 ±0.11 


l b 




M 


219.35 ± 32.90 


1.86 ± 0.15 


69.54 


0.34 ± 0.12 


0.67 


4L38 


L 


13.95 ± 2.09 


2.41 ± 0.20 


65.06 


0.29 ±0.12 


l b 




K 


4.81 ± 0.72 


3.09 ± 0.20 


66.29 


_C 


l b 




J 


12.38 ± 1.86 


4.32 ± 0.28 


64.25 


0.78 ±0.18 


l b 




I 


12.70 ± 1.91 


5.33 ± 0.22 


64.91 


0.85 ± 0.20 


l b 




H4 


32.74 ± 4.91 


6.20 ± 0.24 


61.84 


0.46 ±0.14 


l b 




H3 


38.49 ± 5.77 


6.60 ± 0.21 


60.79 


0.53 ±0.15 


l b 




H2 


35.44 ± 5.32 


7.21 ±0.18 


58.47 


0.84 ±0.19 


l b 




HI 


13.71 ± 2.06 


8.20 ± 0.24 


63.77 


0.88 ± 0.20 


l b 




G 


5.03 ± 0.75 


9.56 ± 0.07 


59.78 


0.52 ±0.14 


I b 




E3 


1.97 ±0.30 


11.42 ±0.28 


66.32 




l b 




El 


4.99 ± 0.75 


15.30 ±0.46 


64.80 


2.10 ± 0.43 


l b 





a The PA is measured from north through east 
k Axis ratio fixed at 1 
c Undetermined by the fit 
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TABLE 3 
Kinematics 



Component 


# of epochs 




i a \ } \ } 


Peak Flux 


Ejection 






[mas yr" 1 ] 




[Jy] 


Epoch 


B 


4 


1.74 ±0.04 


5.7±0.1 




1994.07 ±0.06 


C 


4 


1.4±0.1 


4.6±0.3 




1994.60 ±0.19 


E 


14 


1.00 ±0.02 


3.26 ±0.07 


1.46 


1996. 10 ±0.03 


El 


12 


1.69 ±0.04 


5.5±0.1 


0.30 




E2 


7 


1.29 ±0.06 


4.2±0.2 


0.13 




E3 


12 


1.22 ±0.05 


4.0±0.2 


0.47 




E4 


7 


0.86 ±0.05 


2.8±0.2 


0.12 




F 


4 


0.64 ±0.07 


2.1 ±0.2 


1.98 


1996.13 ±0.16 


G 


9 


1.83 ±0.02 


5.97 ±0.07 


0.24 


2000.54 ±0.03 


H 


4 


1.73 ±0.02 


5.64 ±0.07 


0.43 


2001.43 ±0.02 


HI 


4 


1.9±0.2 


6.2±0.7 


0.03 




H2 


4 


1.71 ±0.11 


5.6±0.4 


0.10 




H3 


4 


1.5 ±0.1 


4.9±0.3 


0.11 




H4 


3 


1.3±0.2 


4.2±0.7 


0.10 




I 


6 


1.3 ±0.1 


4.2±0.4 


0.08 


2001.77 ±0.21 


J 


6 


1.2±0.1 


3.9±0.3 


0.10 


2002.30 ±0.24 


K 


5 


1.1 ±0.1 


3.6±0.7 


0.18 


2003.08 ±0.22 


L 


4 


1.1 ±0.2 


3.6±0.7 


0.11 


2003.74 ±0.22 


M 


3 


1.5 ±0.2 


4.9 ±0.7 


0.39 


2004.50 ±0.12 


N 


3 


0.9 ±0.2 


2.9 ±0.7 


0.54 


2004.57 ±0.19 



